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Abstract. We investigate the chiral phase transition in the quark-meson effective model using optimised perturbation theory 
to one loop. Certain terms in the free energy are frequently omitted in calculations, on the assumption that their contribution 
is negligible. We show that this is not necessarily the case, and that the order of the phase transition, as well as the critical 
temperature, depends heavily on which contributions are included. 
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INTRODUCTION 

The QCD Lagrangian with Nf quark flavours has an 
SU(N/)l x SU(Nf)R chiral symmetry which is sponta- 
neously broken by the formation of a chiral condensate 
at low temperature and baryon density. The location and 
nature of the transition from the phase with broken chi- 
ral symmetry to the chirally restored phase is a topic of 
active research. 

One method of studying the chiral phase transition is 
by using effective models — models which are simpler 
than QCD, yet share QCD's symmetry breaking pattern. 
Examples include the linear sigma model, NJL models, 
and the quark-meson (QM) effective model that we will 
use in the present work. 

Effective models are by their very nature simplifica- 
tions of a more complicated theory. Even so, one often 
makes further simplifications within the models them- 
selves. One may for instance decide to simply omit cer- 
tain parts of the model that don't seem immediately im- 
portant for the problem at hand. As an example, in the 
QM model one often neglects the vacuum contribution 
to the free energy from the quarks, since the chiral sym- 
metry breaking takes place in the meson sector. How- 
ever, one should perhaps not make such assumptions 
light-handedly. It was recently shown that the choice of 
whether or not to include the fermion vacuum contribu- 
tion actually changes the order of the phase transition [ 1 ] . 

In the following, we will look at various approxima- 
tions to the free energy of the QM model, and discuss 
their impact on the phase diagram. 

We will calculate the free energy using optimised per- 
turbation theory (OPT) [2], which is a framework for 
reorganising the perturbative series to obtain improved 
convergence. The present discussion will be limited to 
the one-loop case, as the two-loop calculation is still a 
work in progress. 



THE MODEL 

We consider the quark-meson effective model with 
N scalar fields (01, <p2, . . . , <Pn) and Nf massless quark 
flavours The Euclidean Lagrangian can be written 

as 
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where the three terms describe the mesons, the quarks, 
and the interactions between them, respectively. They are 
defined by: 

1 1 A 

■SO* = -(d^ i )(d^ i ) + -m 2 W i + -(([>4 i ) 2 ,(2) 

^ q = V(Y^-liY4)Vf, (3) 
J% = gWiW^i- (4) 

Here, we have defined T = (1,-/751), where T = 
(ti, Vz, T3) are the Pauli matrices 1 . ji is the quark chemi- 
cal potential, which we assume to be equal for all quark 
flavours, and which is therefore related to the baryon 
chemical potential by ji = jJ.s/3. 

When N = 4 and Nf = 2, this Lagrangian has a 
SU(2)i x SU(2)g symmetry which, assuming m\ < 0, is 
spontaneously broken down to SU(2)v at low energy 2 . In 
the broken phase, (j) acquires a nonzero expectation value 
v, called the chiral condensate, and we therefore write 

= (CT+V,ffl,^,...,7TiV_l). (5) 

In optimised perturbation theory, we add and subtract 
quadratic terms in the Lagrangian: 

& -^^+\x{(y 2 + mm) - l -x{o 2 + mm). (6) 



This only makes sense when N = 4, which is the value we will use in 
numerical calculations. 

2 In the scalar sector, the symmetry manifests itself as 0(N), sponta- 
neously broken down to 0(N — 1). 



The terms with positive sign are treated as mass terms, 
while the terms with negative sign are treated as two- 
particle interactions that give rise to a set of new Feyn- 
man rules. In this way we obtain a systematic resumma- 
tion of selected diagrams from all orders of the naive per- 
turbation expansion. 



FREE ENERGY 

At tree level, we find that the only contribution to the free 
energy comes from the meson sector: 



* = \m 2 v 2 + ^v 4 . 



(7) 



At one loop, after calculating the relevant Feynman 
diagrams in dimensional regularisation and performing 
renormalisation according to the MS scheme, we identify 
the following contributions to the free energy: 

.^i - ^opt + &r + &i + &r (8) 

The first one, J?opt> arises from the OPT interaction 
term, and is defined as 
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Next, we have the vacuum (i.e. T = 0) contribution from 
the scalar fields: 
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Here, A is the MS renormalisation scale. 

The thermal (i.e. temperature-dependent) contribution 
from the meson sector is given by 
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where (O a {p) = y/p 2 + m 2 and ng is the Bose-Einstein 
distribution function. 

Finally, we have the vacuum and thermal contributions 
from the quark sector, which are given by, respectively, 
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dp 



x [n F (co q (p)+lx)+n F (co q (p)-n)} (13) 

where n F is the Fermi-Dirac distribution function and m q 
is a dynamically generated quark mass, m q = gv. 





v=0 


v^O 





200 400 600 800 1000 1200 1400 

Baryon chemical potential [MeV] 
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RESULTS 

The optimisation parameter % is in principle completely 
arbitrary, and if we could carry out calculations to all 
orders, the result would be independent of %■ As it stands, 
we are forced to truncate the loop expansion at a certain 
order, and we must therefore find a prescription for %. 

To that end, we consider the gap equation for the chiral 
condensate v, obtained by differentiating the free energy: 
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Here, IT is the one-loop self-energy IT — 

We now apply a fastest apparent convergence (FAC) 
criterion to %: We choose % so that the perturbative 
correction Fli is as small as possible. In fact, we can 
require that it vanishes completely, TIi = 0. Then it is 
only a matter of solving Eq. (14), now simplified to 
vm 2 = 0, for v. 

We will now look at some common approximations 
to the free energy each obtained by neglecting some 
of the terms in Eq. (8), and we will demonstrate the 
effect of each approximation on the structure of the phase 
diagram. Lattice calculations on two-flavour QCD have 
placed the phase transition at T c ^150 MeV for jj,B = 0, 
and shown that it is of second order. This can, to some 
extent, be used as a measure of the success of the various 
approximations — at least qualitatively. 

Starting with the simplest case, we may say that we 
are only interested in thermal effects, and thus ignore the 
vacuum contributions J?"™ c and j?™ c . Furthermore, it is 



common to use a mean-field approximation to the meson 
sector, in which we ignore the thermal fluctuations, ^ ', 
as well. We show the phase diagram one obtains with this 
approximation in Fig. 1. The solid line indicates a first- 
order phase transition, and the critical temperature for jj,B 
is rather high, T c = 229 MeV. 

For Hb = 0, it was recently shown by Skokov et al that 
if one, in addition to the above, includes the fermion vac- 
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uum term, ^J ac , the order of the phase transition actually 
changes — it becomes a second-order phase transition [1]. 
We have found this to be true for nonzero fis as well, as 
shown in the phase diagram in Fig. 2. The dashed line in- 
dicates a second-order phase transition. The critical tem- 
perature doesn't change much compared to the above. 

Next, we again disregard the fermion vacuum fluctu- 
ations, and instead we include the thermal contibution 
from the meson sector. This will allow us to see the effect 
of taking the mean-field approximation in the first place. 
The phase diagram is shown in Fig. 3, and it agrees very 
well with the results of Bowman and Kapusta [3]. we see 
that the phase transition is first order, and that the critical 
temperature is now significantly lower. Now, T c = 132 
MeV, which is much closer to the QCD critical tempera- 
ture. 

Finally, in Fig. 4, we show the phase diagram when all 
contributions to the free energy are taken into account. 
The phase transition is first order for all jug, and the 
critical temperature is T c = 142 MeV. 



SUMMARY AND OUTLOOK 

We started with a simple approximation, where we only 
included the tree-level contribution from the bosons and 
the thermal contribution from the fermions, and got a 
first-order phase transition with a rather high critical 
temperature. 

We then saw that we could add the vacuum contribu- 
tion from the fermions, which resulted in the phase tran- 
sition becoming a second-order phase transition, or we 
could add the thermal contribution of the bosons, which 
pushed the critical temperature down into the vincinity 
of the QCD critical temperature. 

It would be tempting to guess that by including all con- 
tributions, we would get a second-order phase transition 
with a low critical temperature, but as we have seen, that 
is not the case. When the vacuum fluctuations of both 
fermions and bosons are included, they have a signifi- 
cantly weaker impact on the phase diagram. Part of the 
reason for this is the fact that the vacuum energy of the 
mesons has opposite sign from the vacuum energy of the 
quarks, so to a certain degree they simply cancel out. 

The next step is obviously to include the next loop 
order. At two-loop order we expect to capture more 
aspects of the QCD phase diagram, in particular we 
expect the phase transition to become second order at 
high T and low jug. This is work currently in progress. 
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